Doubly resonant ultrachirped pulses 
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Ultrachirped pulses for which the frequency chirp is of the order of the transition frequency of a 
two-level atom are examined. When the chirp is large enough, the resonance may be crossed twice, 
for positive and negative quadrature frequencies. In this scenario the analytic signal and quadrature 
decompositions of the field into amplitude and phase factors turn out to be quite different. The 
corresponding interaction pictures are strictly equivalent, but only as long as approximations are not 
applied. The domain of validity of the formal rotating wave approximation is dramatically enhanced 
using the analytic signal representation. 

PACS numbers: 32.80.Qk,42.50.-p 



I. INTRODUCTION 

Taking the parameters of a system to extreme values 
beyond their standard domain is a common and fruit- 
ful exercise in physics. It frequently discloses different 
regimes, properties or qualitative changes in the system 
behavior. Appropriate modifications in the theoretical or 
experimental treatments or plainly new tools may also be 
required, and often the new physics found leads to unex- 
pected applications. Examples of much current interest 
in atomic, molecular and optical science are ultracold 
temperatures, or ultrashort and ultrastrong field-matter 
interactions. In this paper we want to put forward and 
examine the limit of ultrachirped pulses. By an "ultra- 
chirped pulse" we mean one where the instantaneous fre- 
quency change is of the order of the transition frequency 
of a two-level (actual or artificial) atom. We shall de- 
scribe in particular linearly chirped pulses that excite the 
resonance twice, at positive and negative instantaneous 
quadrature frequencies. A consequence is a strong diver- 
gence between the analytic-signal and the quadrature- 
model splitting of the field into amplitude and phase fac- 
tors, which causes differences between the corresponding 
formal rotating wave approximations. 

The rotating wave approximation (RWA) is a widely 
applied simplification to treat radiation-matter interac- 
tion neglecting rapidly oscillating terms. Its validity has 
been analyzed for different applications and the de- 
viations from exact results are currently of much interest 
due to the increasing ability to manipulate interaction 
parameter values in different physical settings, and to 
produce strong couplings and/or ultrashort pulses. 

Except for monochromatic, constant-intensity fields, 
the choice of amplitude and phase to describe the field 
is not unique, and it affects the definition of the Rabi 
frequency D,ji{t), instantaneous frequency uj{t), and de- 
tuning A(t). This leads to different interaction pictures 
and different accuracies for the corresponding rotating 
wave approximations (RWA). We demonstrate that the 
amplitude-phase partition provided by the analytic sig- 
nal theory enhances significantly the domain of validity of 



the RWA with respect to the commonly applied quadra- 
ture model partition. To make the paper self-contained 
we shall first review briefly in the following subsections 
the elements of interaction pictures for time dependent 
parameters and of analytic signal theory. Section II is de- 
voted to our study linearly ultrachirped Gaussian 
pulse. 



A. General formulation of intermediate pictures 

Interaction pictures (IP) or "representations" 0, in- 
termediate between the Schrodinger and Heisenberg pic- 
tures, may be most generally formulated in terms of 
a unitary operator Uo{t) 0, which defines, from the 
Schrodinger picture state i^sit), the IP state V'/(i) — 
UQ{t)ips{t)- 'ipii't) evolves according to the dynamical 
equation ihdti^i{t) = Hi{t)^/ji{t), where 



Hiit) = U^it)H'it)Uoit), 
H'{t) = Hit)~Hoit), 
Ho{t) = =ihUoUl 



(1) 
(2) 
(3) 



H{t) is the Schrodinger picture Hamiltonian, and the 
dot denotes derivative with respect to t. In this gen- 
eral formulation the primary object is the unitary oper- 
ator Uo{t), and HQ{t) is "derived" as the Hamiltonian 
for which Uo{t) is an evolution operator. In many text- 
books the emphasis is the opposite: the starting point 
is a splitting of the Hamiltonian H{t) = i?o(0 + H'{t), 
and then Uo{t) is defined as the evolution operator of 
Ho{t) 0. Quite frequently Ho{t) is time-independent 
and Uo{t) = e~*^o*/'', as in typical applications of time- 
dependent perturbation theory, but we stress that this is 
only a particular case and by no means necessary. 



B. Interaction pictures for a two-level atom 

We shall assume a semiclassical description of the in- 
teraction between a laser electric field linearly polarized 
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in i-direction, E(t) = E{t)x, and a 2-level-atoni. Here- 
after we shall refer to E{t) = Eo{t) cos[9{t)] as "the field" . 
In the electric dipole approximation, and for the gen- 
eral case where the transition frequency ujQ(t) may de- 
pend on time, the exact Hamiltonian of the atom in the 
Schrodinger picture is 



H.it) = ^[|2)(2|-|1)(1| 



2 



[(|2)(l| + |l)(2|)(e 



ie{t) 



m 



)],(4) 



This is a "formal" RWA applied blindly at this point, 
without analyzing the frequency content of the neglected 
phase factors, to be distinguished from a more accurate 
treatment in which the phase factor is Fourier analyzed 
and filtered. The term "RWA" refers here to this approx- 
imation and, as we shall see, this crude, formal approach 
will be valid or not depending on the field partition cho- 
sen. 



C. Phase and instantaneous frequency of a signal 



where 9{t) = uj{t')dt' and uj{t) = 9(t) is the instan- 
taneous frequency. Hr = ^R{t) ~ cEa{t) is the Rabi 
frequency with c a real constant, the component of the 
dipole moment in the polarization direction divided by h. 
Note that the definitions of uj{t) and fi/j are not unique, 
as emphasized below. It is usually convenient to write 
the phase as 0{t) = uj^t + if{t), taking as a constant, 
the carrier frequency, so that uj{t) = ujl + (p{t). 

The laser-adapted interaction picture, also called in 
this context rotating frame, based on defining 



Fo(t)=i?L(t) = ^(|2)(2|-|l)(l|) 



(5) 



and 



^-r Hi^{t')dt' /h 

g-^9(t)/2|2)(2|+e*«W/2|^)(l|, 



Uo{t) 

produces the IP Hamiltonian^ 
^A^J - 2 (l + e2*«W)rifl(0 A{t) 



(6) 



(7) 



where we use the vector basis |1) = 
and 

A{t) = Luoit) - uj{t) 



, |2> 



(8) 



is the detuning between the transition frequency and the 
instantaneous frequency. If cjo(0 does not depend on 
time and ujq = ljl, as we shall assume hereafter. 



(9) 



We summarize here some relevant elements of signal 
theory A general real signal can be written in the 
form 



s{t) = a{t) cos e{t), 



(11) 



where a{t) is the amplitude and 6{t) the phase. This 
decomposition of s(t), however, is not unique. In other 
words, we may find different amplitude and phase func- 
tions a' and 9' satisfying s{t) — a'{t) cos6'{t). 

It is customary to define a corresponding complex sig- 
nal. 



Z{t) = A(i)e^®(*) 



with real part s{t), 



Z{t) = s(t) +^s,(^), 



(12) 



(13) 



and imaginary part Si{t) defined in different ways. For 
a given complex signal the instantaneous frequency is 
defined as the derivative of the phase, uj{t) — 9(t), so it 
depends on the imaginary part chosen. 

Moments of the signal may be defined as {t") = 
J s^f'dt/ f s'^dt, and the variance is CTj^ = {t'^)-{t)^. This 
extends to complex signals as well changing — )■ \Z\'^. In 
frequency space, spectral moments are defined similarly 
in terms of the spectrum 
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V2^ J- 



s{t)e-"'^'dt, 



(14) 



and the bandwidth cr^j is defined as the root of the vari- 
ance. For a real signal 5'(— a;) = '5'*(i^) and [^(ti;)^ is 
symmetric about the origin, so (w) = 0. It is also possi- 
ble to extend the spectrum concept to the complex signal, 
or to parts of it, and break this symmetry. 



Applying now a RWA to get rid of the counter-rotating 
terms we end up with 

= 2 17;^(t) A{t) )■ (10) 



This was termed, with a different notation "quasi-interaction pic- 
ture" in In fact, according to the definition of IP given above, 
it is a perfectly canonical one. 



1. The quadrature signal 

From the form (|lip . it is natural to write the complex 
signal as 

5,(0 =a(t)e'«(*), (15) 

which is called the quadrature (model) signal. As Eq. 
(fTTj) is not unique, the quadrature signal is not unique 
either. Sq{uj) is the corresponding spectrum. 
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2. The analytic signal 



The analytic signal is a peculiar complex signal chosen 



(16) 



The imaginary part Sa,i{t) is the Hilbert transform of the 
real signal s{t). Sa{t) can also be written in terms of its 
amplitude Aa{t) — (s^ + J^^^ and its polar phase 9ait) 
as 



(17) 



One of the advantages of the analytic signal is that it 
puts the low frequencies in the amplitude and the high 
frequencies in the phase factor e**"^*-* Q. If the spec- 
trum of the amplitude is of finite support, the support 
of the spectrum of the phase factor does not overlap and 
this makes the pair [Aa{t),9a{t)] unique. The spectrum 
of Sa{t) is Sa{uj) = Sq{uj) + S*{-uj) if w > and zero 
otherwise. 

The maximum possible deviation of a quadrature sig- 
nal from the analytic signal, at time t, is given by 



\Sa{t) - Sg{t)\ < 



\Sq{LLl)\duJ. 



Another criterion of closeness is that 



\Sa{t)-Sg{t)\''dt^2 \Sg{Lj)\^dOJ 



(18) 



(19) 



should be small, i.e., the spectrum of the quadrature sig- 
nal should be predominantly in the positive frequency 
domain [sj. 



II. RWA WITH THE QUADRATURE AND THE 
ANALYTIC SIGNALS FOR A LINEARLY 
CHIRPED GAUSSIAN PULSE 

We consider now an electric field with a Gaussian en- 
velope and a linear chirp. 



(20) 



as the real signal, where i?o(0 — foe^°* is its ampli- 
tude. Not to restrict the numerical examples to a par- 
ticular system or region of the spectrum, we shall from 
now on use dimensionless variables. We define i = a^^^t, 
Cb = Lo/a^/"^ , and apply the same scaling to all times and 
frequencies. Similarly, b = b/a. From the definition of 
the dimensionless Rabi frequency, = Qn/a^^'^, the di- 
mensionless amplitude of the field is Eq — Eqc^* , where 
So — c£o/a^^'^. Dimensionless Hamiltonians are defined 
as H = H / {a}/'^h). To avoid a heavy notation we finally 



drop all the tildes hereafter. The dimensionless version 
of Eq. (HOI) is thus 



E{t) = foe"*' cos [ujLt + bt^), 



(21) 



in terms of a carrier or central frequency wl, a chirp 
parameter 6, and maximum amplitude E^. In these units, 
all pulses have equal duration, with at — 1/2. 




FIG. 1: Linearly chirped Gaussian pulse. Parameters: b = 2, 
£q = 4^, and cjl = 10^. 

The corresponding quadrature model signal is 

Eq{t) ^ £oe-'" e'^'^^'+^'"\ (22) 

with amplitude Eoq{t) = foe~*^, phase 9q{t) = ujLt + bt"^, 
and instantaneous frequency (^q{t) — ujl + '2bt corre- 
sponding to a linear chirp. The dimensionless interaction 
picture Hamiltonian in the rotating wave approximation 
(fTOl) becomes 



E[l,RWA,qit) 



-Aq{t) nRq{t) 
riliqit) Aq{t) 



(23) 



where Aq{t) = —2bt if ujQ{t) — ojl. 

The analytic model signal provides instead the com- 
plex field Ea (t) according to Eq. ([Tl 



Eait) = Eoait)e 



(24) 



The corresponding interaction picture is different from 
the one using the quadrature model, and the RWA Hamil- 
tonian pH)) becomes 



Hl,RWA.a{t) = 



1 



ilRait) Aa{t) 



(25) 



where Aa{t) = u}o{t) - ^^ait), with uiait) = 6a{t 
To calculate Ea{t) we may use the relation 

Eiuj) ^ ^[Eq{u;) + E;{~u;)], 



(26) 



where E{uj) is the spectrum of E{t), and 

1 Sge^i'^~''^Ly/''^(^-i'>) 



Eqiiu) 



/2tt J- 



Eq{t)e 



V2(i - 



(27) 
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Then, from Eq. ([16]), 
So 



Ea{t) 



where z = z{t) = ty/1 — ib — iw_L/(2v'l — ib) and w{z) = 
erfc{—iz) is the Faddeyeva or w-function Now 
the amplitude and the phase are defined as EQa{t) = 

^E^{t) + {lm[Ea{tW and Bait) = Im[ln£;,(i)]: respec- 
tively. 

From Eq. (HH) [|, and for our Gaussian electric field, 
Eq{t) ^ Ea{t), at least in the central part of the pulse, 
for 



(29) 



where aujq = \/ (f ^ + is the bandwidth of Eqit). 

4i 




FIG. 2: (Color online) (a) Phase of the analytic signal 9a{t) 
from — TT to TT. (b) The instantaneous frequency uJa{t) (solid 
red line) following the linear form of tOqit) = uiL-\-'2bt (blue line 
with dots) in the central interval and tending asymptotically 
to zero. Parameters: b — 2, So = 4-\/2, and ljl = f0\/2- 



pictures, although different, are all equivalent since they 
lead to the same Schrodinger dynamics and to the same 
populations. Applying the RWA gives, however, different 
results. 

Considering that the atom is initially in the ground 
state, we have solved numerically the system of differen- 
tial equations for the wave function amplitudes with the 
exact Hamiltonian ([T]), and approximate ones (1231) and 
(113. 

The validity of the RWA is usually discussed for time- 
independent frequencies and linked to assuming first the 
weak-coupling condition [1, Q , 



(30) 



and then a quasi- resonance condition jAj/wo ^ 1- This 
later condition is clearly violated out of resonance for ul- 
trachirped pulses. Nevertheless, if the condition (pO| is 
satisfied and ujl is also large enough to satisfy the close- 
ness condition (|29p . the populations P2{t) driven by the 
two approximate Hamiltonians are equal to each other 
and to the population driven by the exact Hamiltonian. 
When ([m is not satisfied, however, the populations P2{t) 
driven by these approximate Hamiltonians are different 
and the range of validity of the RWA is much wider for 
the analytic signal than for the quadrature, as we shall 
see. 



1 

0.8 
0.6 



0.4 

0.21 





-0.4 -0.2 




t 



0.2 0.4 



We have defined 9q{t) — uj^t + bt^ to have a linear 
chirp, but 9ait) will have in principle a different form. 
When Ea{t) ^ Eq{t) in a central interval, 9a{t) follows 
the quadratic form of 9q{t) there, but it becomes constant 
well before and after, see Fig. [2][a). The "parabola" is in 
fact cut into pieces corresponding to the principal branch 
of the logarithm, from — tt to tt. Similarly, Fig. [2{b) 
shows that the instantaneous frequency ujait) follows the 
linear form of Wg(t) = uj^ + 26t in the central interval 
but it tends asymptotically to zero. The details of the 
transition are discussed in the appendix [X] 

A. Population inversion 

If the dynamics is performed exactly, different ampli- 
tude and phase splittings of the field produce the same 
results. Only the real field matters, and the interaction 



FIG. 3: (Color online) Population of the excited state for 
the exact dynamics and for the RWA approximations with 
the quadrature signal and the analytic signal. The carrier 
frequency is large enough so that the three lines essentially 
coincide. Parameters: 6 = 3500, So = 200, and cjl = 20000. 

The parameters of the system, considering that wo(t) = 
wl, are £q, ujl, and b. For the values So = 200, b = 3500 
and ujl — 20000, see Fig. [3l the population is fully trans- 
ferred to the excited state according to the exact dynam- 
ics. We then decrease ujl up to 700. Figures |3]l5] show 
that in this range of carrier frequencies, the population of 
the excited state P2 (t) with the analytic signal RWA fol- 
lows the exact results, while the quadrature RWA fails. 
At UJL — 700, see Fig. El the quadrature RWA model 
predicts full inversion, whereas the analytic RWA model 
and exact results give a complete population return. 

This behavior shows up that the condition (PH)) for 
the quadrature RWA is much more strict than for the 
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FIG. 4: (Color online) Population of the excited state for 
(a) the exact dynamics (b) the Hamiltonian in the RWA for 
the quadrature signal (blue line with dots) and the analytic 
signal (solid red line). Compare to the previous figure. As 
ljJl decreases, only the analytic signal RWA follows the exact 
dynamics. Parameters: b — 3500, £o = 200, and ujl = 5000. 



analytic signal RWA. To reach the agreement between 
the quadrature RWA and the exact dynamics we have to 
set ujl = 20000, 100 times greater than the maximum of 
V,iiq{t), r2flq(0) — 200, whereas the analytic signal RWA 
still follows the exact results when wl = 700, only 3.5 
times greater than the maximum, Jlija(O) = 200.08. For 
~ 600 (figure not shown), a small discrepancy appears 
between the populations of the exact and analytic RWA 
model which increases with decreasing ujl- 

The complete population return shown in Fig. [S] is 
stable, for b between b — 2000 to 6 = 4000. Since vary- 
ing b amounts to shifting the NF resonance, this provides 
"windows" of excitation of controllable duration, poten- 
tially useful for optical gating, with switching times much 
shorter than the pulse duration. 

There is, according to the construction of the pulse 
and the condition ujl = '^Oi a resonance region at t — 0. 
The early excitation at negative times in Figs. |3][5]is due 
to the crossing of an additional resonance region where 
the instantaneous frequency becomes, in the quadrature 
model, —Wo, see Fig. [51 Using as before wq = wl, we have 
that \u!l\ = <jJl + 26t has solutions i = 0, corresponding 
to the nominal resonance of the pulse, and also 



t = tr —ujL/b. 



(31) 



As shown in Figs. |4][5l the quadrature RWA does not no- 
tice this "negative-frequency" (NF) resonance, since the 
negative instantaneous frequency corresponds in fact to 
a large quadrature detuning Ag ~ 2ajo at tr- Instead, the 
instantaneous frequency for the analytic signal is positive 
at tr and gives = 0, see Fig. [S] The condition for the 
NF resonance to occur out of the pulse (therefore having 
no consequence) is that tr » 3crt, i.e.. 
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FIG. 5: (Color online) Population of the excited state for: (a) 
exact dynamics; (b) RWA for the quadrature signal (blue line 
with dots) and the analytic signal (solid red line). Parameters: 
b = 3500, £o = 200, = 700. 
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UJL » 3&/2. 



(32) 



FIG. 6: (Color online) Instantaneous frequencies ujq{t) = 
UJL + 2bt (blue line with dots) and uja{t) (solid red line). Pa- 
rameters: b = 3500, So = 200, and ujl = 5000. The dashed 
line is at ujl- The arrows indicate the early, and nominal 
resonance regions for the pulse, when uja = ujl- 



This is independent of the magnitude of the interaction. 
Since the difference between the single-resonance condi- 
tion ([5^ and the closeness condition is just a numer- 
ical factor, when 6 >> 1, the separation of analytic and 
quadrature approximate RWA dynamics is associated 
with the occurrence of this early NF resonance within 
the pulse. Fig. [7] depicts an example of an ultrachirped 
real signal capable of inducing a double resonance. The 
passage through zero instantaneous frequency is evident 
as a gap in the lower envelop and nearby slow oscillations. 



B. Adiabatic approximation 

It is interesting to complement the above results by ex- 
amining the validity of the adiabatic approximation. The 
condition for the states evolving with and to 
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FIG. 7: Ultrachirped signal. The arrows indicate the reso- 
nance regions where u}l = ^a- Parameters: 6 = 30, £o = 5, 
and u}L =40. 



alization of the parameters of Fig. [5] could be as fol- 
lows: = 27r • 1.019 GHz, a = {2iTf • 0,21 • lO^^ Hz^; 
h = (27r)2 . 7.4 • lO^'"^ Hz^, which gives at = 54.652 ns, and 
a^q = 2tt ■ 5.09 GHz. 

While this work has been essentially a curiosity-driven 
exploration, applications may be envisioned in state de- 
termination, optical gating and interferometry, since the 
timing of the resonance regions can be controlled. Spe- 
cific fields we consider for future work are: Circuit Quan- 
tum Electrodynamics, in which the two-level system and 
its interactions are highly tunable and potentially time 
dependent, and also cold atoms in counterpropagating 
laser beams inducing a Raman transition |12| . 



behave adiabatically following the instantaneous eigen- 

states is 



\nd(t)\ » \\nAd{t)l 



(33) 



where the subscript d = q,a refers to the amplitude-phase 
decomposition, quadrature or analytic signal, Vtd{t) = 
V^fid + effective Rabi frequency, the instan- 



taneous energies for the Hamiltonians ([25]) and ([25]) 

are Ed±{t) = Vtd{t)/2, and ilAd{t) = [^Rd^d{t) ~ 

/:^d{mRd{t)]/nl{t). 

In our numerical example, according to Eq. (1291) . 
lol = 20000 is large enough to make the analytic and 
the quadrature signal models very similar. For this ujl 
the adiabaticity condition is satisfied both for Eqit) and 
for Ea{t). As we diminish uj^, however, the adiabaticity 
condition is still satisfied for Eg, but not for Ea{t). This 
is quite evident e.g. in Fig. |4l note the smooth curve 
of the quadrature signal model versus the coherent tran- 
sients of the exact results or the analytic signal model 

mm- 



III. DISCUSSION. CONCLUSIONS 

In this paper we have explored the consequences of 
going beyond common excitation regimes of a two-level 
quantum system, in particular we have seen that for ul- 
trachirped fields the resonance condition may be satisfied 
twice, for positive and negative quadrature instantaneous 
frequencies along the linear chirp. The formal rotating 
wave approximation is more robust by using analytic sig- 
nal theory for the complex signal and the corresponding 
interaction picture. The reason is that the analytic signal 
reinterprets the chirp, assigning positive frequencies to 
the pulse region with negative frequencies in the quadra- 
ture model. 

Numerical examples demonstrate that the necessary 
dimensional parameters to see double resonances are 
within reach in the microwave domain with the cur- 
rent technology of pulse generators. A dimensional re- 
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Appendix A: Asymptotic values of the phase of the 
analytic signal 

For a linearly chirped Gaussian pulse the phase 9a {t) — 
lm[\nEa{t)] of the analytic signal ([251) tends to a con- 
stant, see Fig. [2| 2; = z{t) = Wl — ib — iuJL/i^V^ — ib) 
is a linear function of t, so as \t\ increases z{t) becomes 
larger in modulus and asymptotic expressions of the w 
may be used. The asymptotic behavior of Ea {t) for large 
\z\ and Im[z(t)] > is 



-u:l/4{l+zb) 



+ 2e""'/''(^+*'')e~(^')'(*) 
whereas for Im[z(t)] < [isj . 



(Al) 
(A2) 



/TTZ*{t) 



(A3) 



Then, 



Bait) ^ arctanj ^ ^ [,^(^)f | + (A4) 

where the negative and positive signs correspond to 
Im(2;) > and Im(z) < 0, respectively, n = 0,1,2,... 
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FIG. 8: (Color online) In {\Ea{t)\) (red triangles) is essentially 
equal to ln(|(72(t)|) (dashed blue line) in the central interval, 
and to In (|(7i(i)|) (green solid line) in the asymptotic regions. 
Parameters: b — 2, £o ~ 4^2, and ujl — 10-\/2. 



Therefore, the asymptotic expansion for the phase of the 
analytic signal will be 



2?™. 



(A7) 



For n = 0, the phase remains inside the principal branch 
of the hiEa{t), see Fig. [21 As the phase reaches these 
constant values, uja{t) when t ±oo. 

To identify the transition between the central pulse 
and asymptotic regimes let us compare different terms 
with the full expression. Defining 



51 W 



So 
2 



le 



-c^i/4(l~»6) 



-"!/4(l+»fc) 



(A8) 



and 



with 



2(1 + 62) 1/4 cos a + ibt+^) sin a] 

^-c^i/4(l+&2) 



bujl 1 

' - arctan o. 



4(1 + 52) 2 



(A5) 



(A6) 



Independently of the signs ± in Eq. (jA4p . when t ±oo, 
the argument of the arctangent tends to +oo or to — oo. 



52 (i) 



_ r £oe-'^'/4(i+'f')e-(^*)'(*) Im[z{t)] > 



foe 



Im[z{t)] < 



(A9) 



and comparing In [|i?a(t)|] with In [|.gi(i)|] and In [|52(0I]j 
we see that (jA9p dominates during the central interval, 
and (lASp in the outer time regions, see Fig. [51 Transi- 
tion regions with interferences occur near the times when 
ln[|gi(t)|] and ln[|g2(t)|] are equal. For the parameters 
b = 2, So ^ AV2, and lul = 10\/2, see Fig. [HI these 
instants are ti = —3.862 and t2 = 3.596, so the pulse is 
essentially within the central interval. 
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